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Abstract 

A family of one-dimensional multi-species reaction-diffusion processes 
on a lattice is introduced. It is shown that these processes are exactly 
solvable, provided a nonspectral matrix equation is satisfied. Some general 
remarks on the solutions to this equation, and some special solutions are 
given. The large-time behavior of the conditional probabilities of such 
systems are also investigated. 



PACS numbers: 82.20.Mj, 02.50.Ga, 05.40.-a 



1 Introduction 



In recent years, the asymmetric exclusion process and the problems related to 
it, including for example the bipolymerization jlj], dynamical models of interface 
growth || , traffic models || , the noisy Burgers equation Q , and the study of 
shocks H ||, have been extensively investigated. The dynamical properties of 
this model have been studied in [6-8]. As the results obtained by approaches 
like mean field are not reliable in one dimension, it is useful to introduce solvable 
models and analytic methods to extract exact physical results. Among these 
methods is the coordinate Bethe-ansatz, which was used in || to solve the 
asymmetric simple exclusion process on a one-dimensional lattice. In [jH)|, a 
similar technique was used to solve the drop-push model , and a generalized 
one-parameter model interpolating between the asymmetric simple exclusion 
model and the drop-push model. In |JL2], this family was further generalized 
to a family of processes with arbitrary left- and right- diffusion rates. All of 
these models were lattice models. Finally, the behaviour of latter model on 
continuum was investigated in [T^ | . The continuum models of this kind are also 
investigated in (l4|, In |l6| a generalization of such processes was studied 
which contained annihilation of particles as well. 

In all of these, people have been mainly concerned with the so-called single- 
species processes, in which only one kind of particles exist and move on the 
lattice (or the continuum). Another interesting problem is the study of multi- 
species systems in them several kinds of particles move and interact on a lattice. 
In p7j |, single-species systems have been characterized for them the equations 
governing the evolution of the TV-point functions contain N- or less- point func- 
tions. This has been done for multi-species systems in |lq|. In p9[ , two-species 
reaction-diffusion systems have been introduced that are solvable in the sense 
that the S matrix corresponding to them is factorizable into the two-particle 5* 
matrices. It is found there that the the criterion for this is that the interactions 
(which are of the nearest-neighbor type) are to be so that the S matrix satisfy 
a kind of spectral Yang-Baxter equation. 

We follow the same line. That is, we investigate interactions which can be 
written as boundary conditions for the probability functions. These interactions 
preserve the total number of particles, so that if one begins with N particles, 
knowing the iV-point probabilities is enough to know everything about the sys- 
tem. Using the coordinate Bethe-ansatz, it is found that for this ansatz to be 
consistent, the S matrix should satisfy a kind of spectral Yang-Baxter equation 
[ [l9[ . However, the S matrix is of a special form containing the boundary con- 
ditions (or interactions), and not every solutions of the spectral Yang-Baxter 
equation can be used to construct such a solvable model. We investigate the 
spectral equation the S matrix should satisfy and show that this is equivalent 
to a non-spectral equation for the boundary conditions. This is independent of 
the number of species. 

The scheme of the paper is the following. In section 2, the it is shown that the 
prescription of investigating multi-species reaction-diffusion systems in terms 
of diffusion systems equipped with suitable boundary conditions is studied. In 
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section 3, the Bethe-ansatz solution for such (solvable) systems is obtained and 
its large-time behavior is investigated. In section 4, the solvability criterion is 
obtained and it is shown that this criterion is a nonspectral matrix equation. In 
section 5, some general properties of the solutions of the solvability criterion is 
studied. Finally, in section 6 some special solutions of the solvability equation 
are studied. 

2 Mult i— species reaction— diffusion systems and 
the boundary conditions 

Consider a system consisting of N particles on a lattice, drifting to the right 
with unit rate if the right neighboring site is empty, and interacting with each 
other only if two of them are adjacent. Suppose that there are n kinds (or 
species) of particles and the interaction between the particles is just in the form 
that if two particles A a and A@ are adjacent to each other, they may change to 
A 1 and A s with the rate That is, the allowed processes are 

A a % 0A Q , with rate 1 

A a A p -> A^A 5 , with rate 6$ (1) 

These processes result in the master equation 

P Q1 '- < Q ~ (xt, ■ ■ ■ , x N ; t) =P a ^- > a » (:n - 1, ■ • • , x N ; 1) + ■ ■ ■ 

+ P ai >-> 0W (si,--- ,x N -l;t) 
-NP a ""> a »(x 1 ,--- ,x N ;t), (2) 

if Xi < x i+ i — 1. The symbol P ai '"' ' aN (x\, ■ ■ ■ ,XN]t) denotes the probability 
of finding a particle of type a\ in xi, a particle of type at2 in x 2 , ■ ■ ■ at the time 
t. The so-called physical region consists of the points satisfying Xi < Xi+\. If 
Xi = x i+ i — 1, the interactions change the equation. For clarity, let's write the 
evolution equation for the two-particle sector: 

P af3 {x, x + l) =P af3 {x - 1, x + 1) 

+ E b$F"{x,x + l) 

- P afi {x,x + l), 

where 

B#:= J2 6j (4) 

(lSma/3) 

Defining the diagonal elements of b as 

&3 :=1 - E b l% (5) 

( 7 <5)#(a/3) 



~B a(3 P al3 {x,x + l) 

(3) 
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it is seen that ([|) can be written as 

P al3 (x, x + 1) = P a(3 {x - 1, x + 1) + b^P~< 5 {x, x + 1)- 2P af3 {x, x + 1), (6) 

where summation is implied on repeated indices. Comparing this with (Q). 
it is seen that it can be written as (|J) provided one introduces the boundary 
condition 

P<*P(x,x)=b$Pi s (x,x + l), (7) 

or, in a more compact form, 

\P(x,x))=b\P(x,x + l)). (8) 

The matrix b should satisfy two criteria. First, its non-diagonal elements should 
be nonnegative (since they are rates). Second, the sum of the elements of each 
of its columns should be one. This can be written in a compact form as 

(«|®(*|6=(a|®(*|, (9) 

where 

s a := 1. (10) 

Note that if the number of species is one, the asymmetric simple exclusion 
process |J is obtained. 

3 The Bethe— ansatz solution 

As in fl9fl , one can write a Bethe-ansatz solution for (|J) with the boundary 
condition 

\P(- ■■ ,Xk= x,x k+ i = x, ■ ■ ■ )) = b k . k+ i\P(- ■■ ,Xk = x,x k+ i = x + 1, •••)), 

(11) 

where 

b kM+1 := 1 ® •■ • <g> 1 ® b ®1® ■■ • <g> 1. (12) 

k.k+l 

We take the ansatz 

|P(x;<))=e £t |vI/(x)), (13) 

and it is seen that |^(x)) should satisfy 

E\^{x x ,-- - ,x N )) =\#(xi - l,-- - ,arjv))+--- 

+ \*(xi, ■ ■ ■ , x N - 1)) - N\*{x u ■■■ , x N )), (14) 
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and 

|*(- ■ • ,Xk = x,x k+1 =x, •■■)} = &fc,fc+i|*(- •• ,x k = x,x k+l =x+l, ■■■)). 

(15) 

The Bethe-ansatz is that, one takes the following form for ^(x)}. 

|*(x))=^^e^P)- x |V), (16) 



where \tp) is an arbitrary vector and the summation runs over the elements of 
the permutation group. Plugging this in ( fTi]) results in 



jv 



k=l 



The boundary condition (|15|) yields 



1 _ e-^+^&u+il A a +\l- e^b k , k +i 



A<j<j k — 0. 



(17) 



(18) 



Here a is that element of the permutation group which only interchanges p k and 
p k +\. From this, one obtains 



Aaa k = S ktk+ i[a{p k ),a(p k+ i)]A a , 
where the matrix S is defined through 

S(pi,p 2 ) :=-(!- z x b)- 1 (l-z 2 6), 



(19) 



(20) 



and the definition of S k ,k+i is similar to that of bk, k +i in ([121). we have also 
used the definition 



(21) 



This shows that one can construct A a 's from A\ by writing a as a product of 
tTfc's. But these elements of the permutation group satisfy 



This means that 



CfcCfc+iCfc = ak+icr k a k +i- 



A = A 



(22) 
(23) 



or 



Sk,k+l {Pk+1 , Pk+2)Sk+l,k+2 (Pk , Pk+2 )Sk,k+l{Pk, Pk+l) 
z Sk+l,k+2(pk,Pk+l)Sk,k+l(Pk,Pk+2)Sk+l,k+2{pk+l,Pk+2)- (24) 
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This can be written as 



[S(p 2 ,p 3 ) <B> 1][1 <8> 5(pi,p 3 )][S(pi,P2) ® 1] =[1 <8> S(pi,Pa)][5(pi,P3) <8> 1] 

x[l®S(pa,j9a)]. (25) 

Writing the S 1 matrix as the the product of the permutation matrix II and an 
R matrix: 

Sk,k+i ='■ Rk,k+iRk,k+ii (26) 

5j) is transformed to 

R23(P2,P3)Rl3(Pl,P3)Rl2(Pl,P2) = R\2 (Pi , P2)Rl3 (Pi , ^3^23 (P2, P3) ■ ( 2 7) 

This is the spectral Yang-Baxter equation. 

Provided this condition is satisfied, it is easy to see that the conditional 
probability (the propagator) is 

C/(x; i|y; 0) = J ^pe"*™ £ ^e^W'V E ^ , (28) 

where the integration region for each pi is from [0, 2ir], and we have taken A e = 1. 
(e is the identity of the permutation group.) Note that (0), and the condition of 
nonnegativity of the nondiagonal elements of b ensure that the absolute values 
of the eigenvalues of b don't exceed 1. So there is no singularity in S(pi,f>2) 
except at p\ — 0, and this is removed by setting pj — > pj + it, where one should 
consider the limit e — > + . This is the same as what has been done in || and 
[ [To) , for example. Using this propagator, one can of course write the probability 
at the time t in terms of the initial value of the probability: 

|P(x;t))=^(x;t|y;0)|P(y;0)) (29) 
y 

For the two-particle sector, it is not difficult to obtain U. In fact, as there is 
only one matrix (6) in the expression for U, one can treat it as a c-number and 
the problem is reduced to that of p6[ , with A replaced by b. So, 

U(x;t\y:Q)=e~ 2t 



(xi - yi)! (x 2 - 2/2)! 

e 7T^ T77^ (-!+ • ( 3 °) 

^ (/ + x 2 - yi)! (xi - 7/2)' V xi— 2/2 + I' 



i=0 



One can decompose the vector space on which & acts into a subspace on which 
6=1 (eigenspace of b corresponding to eigenvalue 1) and another invariant 
subspace. This is done by decomposing the unit matrix into two projectors: 

1 = Q + R, (31) 
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where Q and R are projections satisfying 



QR = RQ = 0. 



(32) 



Q is projection on the eigenspace of b corresponding to the eigenvalue 1, and R 
is projection on the other invariant subspace of b. Using this, one can write U 
as 



J7(x;t|y;0) 



-21 



(xi - yx)\ (x2 - V2)\ 

fl+x-z—yx ifi—Vi 



2=0 



(l + x 2 - vi)\ (xi - y 2 )! 



-1 



5Cl - J/2 + 1 



+ 



+e" 



pi— i/i 



E 



(xi - j/i)! (x 2 - y 2 )! 



Here we have used 



(Z + x 2 -yi)\ (xi -y 2 )\ 



b = b(Q + R) = Q + bR. 



t b 



Xl -2/2 + 1 



i?. 

(33) 
(34) 



As the eigenvalues of b, other than 1, are assumed to have moduli less than one, 
the second term in (|3^) is the same as (33) in that is, a term obtained from 
the boundary condition corresponding to annihilation (A < 1 in ]l6| ). The first 
term corresponds to an asymmetric simple exclusion process M . The large-time 
behavior of these two terms are also simply obtained. The large-time behavior 
of the first was obtained in [|l3], and that of the second in |L6|. At large times, 
the second term is found to be independent of b (or A) and vanishing faster than 
1/t. Also, the summation of this term vanishes as t tends to infinity. In fact, 
using p6[ it is seen that 



the second term = -; e 

2?rf 



[(xi- H i-i) 2 +(x2-»2-i) 2 ]/(2t) 
[(x 1 ~y 1 -t) 2 + (x 2 -y2-t) 2 ]/(2t)^ ^ t ^ OO. 



(35) 



So at large times only the first term of ( |33| ) survives. This means that at large 
times, the propagator is proportional to the projection on the eigenspace of b 
corresponding to the eigenvalue 1 (the projection on the equilibrium subspace of 
b) and the proportionality constant is simply the propagator of the asymmetric 
simple exclusion process. 

To conclude, for large times the two-particle conditional probability is that 
of an asymmetric simple exclusion process projected on the eigenspace of b 
corresponding to its unit eigenvalue. 
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4 Solvability criteria for the boundary condi- 
tions 

From (^p|), it is seen that S(px,p 2 ) is a binomial of degree one with respect 
to z 2 '■= e lp2 . Putting this in (^5|), one arrives at a quadratic expression with 
respect to Z3. The coefficients of this expression are, of course, matrices depend- 
ing on z\ and z 2 . It is easy to find the roots of this expression for Z3. In fact, 



putting Z3 = z\ in (25), one arrives at the identity 



[S( P2 , Pl ) ® l][5(pi,pa) ® 1] = [1 <8> S(pi s pa)][l ® S( P2 ,Pi)}. (36) 

(We note that S(pi,p2)S(p2, P i) = 1.) Also, putting Z3 = Z2, another identity 
is obtained: 

[1 ® S(pi,P2)][5(pi,Pa) ® 1] = [1 ® S(pi s P2)][5(pi,Pa) <8 !]• (37) 

These two identities show that the roots of the quadratic expression for Z3 are 
z\ and z 2 . That is, one can write that expression as 

(23 - zi)(z3 - z 2 )Q(zi,z 2 ) = 0. (38) 
So, ( p5| ) is equivalent to Q = 0, which itself is obtained by putting 2:3 = in 

[(1 - Z2&)- 1 ® 1][1 ® (1 - Zl&r X ][(l - ^l6) _1 (l - 2 2 6) ® 1] 

= [1 ® (1 - »i6)- x (l - z 2 b)][{l - Zib)' 1 ® 1][1 ® (1 - z 2 b)- 1 ]. (39) 

Inverting both sides, one arrives at 

[(1 - z 2 b)- 1 (l - zib) ® 1][1 ® (1 - zi6)][(l - z 2 b) ® 1] 
= [1 ® (1 - z 2 b)][{l - zib) ® 1][1 ® (1 - z 2 b)- 1 (l - Zl b)]. (40) 

This is a quadratic expression in terms of Z\. For Z\ = 0, (^) gives the identity 

[(1 - z 2 b)~ l ® 1][(1 - z 2 b) ® 1] = [1 ® (1 - z 2 b)][l ® (1 - Z2&)- 1 ], (41) 

while for z\ = z 2l the identity 

[1 ® (1 - Zib)] [(1 - z 2 b) ® 1] = [1 ® (1 - 226)] [(1 - z 2 b) ® 1] (42) 

is obtained. So, the quadratic expression corresponding to (^) is equivalent to 

zi{z x -z 2 )Q{z 2 ) = 0, (43) 

and to find Q, one simply uses the coefficient of z\ in (flof). This is 

(1 - z 2 bi 2 )~ 1 b 12 b 23 (l - z 2 b 12 ) = (1 - 2:2^23)^12^23(1 - Z2623)" 1 , (44) 
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or 

612623(1 - z 2 b 12 )(l - z 2 b 23 ) = (1 - z 2 b 12 )(l - z 2 b 23 )b 12 b 23 . (45) 

This is a quadratic expression in z 2 . But the coefficients of z\ an d z | are 
identities. So the only remaining equation is 

612623(612 + b 23 ) = (612 + 623)^12623, ( 46 ) 

or 

6i2[6i2,6 2 3] = [6i2,6 23 ]623- (47) 

( [47l ) is equivalent to (|2^). But it is seen that (^) is non-spectral, whereas ( p5| ) 
is spectral. So it is far simpler to seek the solutions to ([47j ) than to seek those 
of©- 

To summarize, a matrix 6, or the reactions ([!]), correspond to an exactly 
solvable reaction diffusion system on a one-dimensional lattice, provided b sat- 
isfies ( f47| ) and (0) (or (||), equivalently), and the non-diagonal elements of b are 
nonnegative. 

5 General properties of the solutions to the solv- 
ability criteria 

Solutions to ( f47| ) enjoy two general properties. First, if b is a solution, then 

b' := ah + 0, (48) 
is another solution for constant a and (3. If b satisfies ^, then 

(s|(g)(s|6' = (a + /3)(s|(g)(s|. (49) 

So, putting (3 := 1 — a, ensures that b' satisfies (^|). If a > 0, then the nondiag- 
onal elements of b' are nonnegative provided the nondiagonal elements of b are 
nonnegative. So 

b' := ab + (1 - a) (50) 

corresponds to a solvable system (for a > 0) if b does. It is easy to see that the 
meaning of this transformation is simply to multiply the reaction rates by a. 
Second, if 6 is a solution to ( f47j ) , then 

b' := u (g) u b u^ 1 (g> u^ 1 (51) 

satisfies (^7|) as well. Here u is an arbitrary (nonsingular) matrix. This trans- 
formation, however, does not necessarily respects the conditions (|^) and non- 
negativity of the rates. So, another problem arises. Suppose b is a solution to 
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(p7|), and we want to obtain a solvable system using the transformation (|T| 
We must have 



(s\ ® (s\u<g>ub= (s| ® (s|w(g)u. (52) 

This means that m must change (s| to some (s'\ so that (s'\ <E> (s'\ is a left 
eigenvector of b corresponding to a unit eigenvalue. One may search in the 
eigenvectors of b to find whether there is an eigenvector of the form (s'\ ® (s'\. 
If there is such an eigenvector, then any matrix u which changes (si to (s'\ can 
be used to obtain b' according to (|5l"l). This b' satisfies (^7|) and (||). But its 
diagonal elements may be nonnegative or not; this should be checked separately. 
If non of the eigenvalues of b are of the form (s' | <8> (s'\, then this method cannot 
be used to obtain a solvable system. This method resembles very much to that 
used in pof . 



6 some special cases 

case I: b 2 = a + 13b (a and (3 are numbers). In this case, one can define 

b':=b + 7 , (53) 

with 7 satisfying 

7 2 + 07 - a = 0, (54) 

to obtain 

b' 2 = {13 + 2i)b'. (55) 
Putting this b' in (^7|), one obtains the braid equation for b': 

& 'l2 & 23 fo 12 = &23 b 12 b 23- ( 56 ) 



From (p5[), it is seen that b' either can be scaled to a projection (b a = b'), or is 
nilpotent. One concludes then that any nilpotent or projection solution to the 
(nonspectral) braid equation is a solution to ([f7|). One can then use any linear 
combination of this solution with the unit matrix as another solution to that 
equation. Note, however, that these solutions of ( ft7| ) do not necessarily satisfy 
other criteria of the solvable system, that is nonnegativity of the nondiagonal 
elements and (||). An inspection of the solutions obtained in |h| shows that 
solutions 1 — 15, and 17 are of this type. As mentioned in the previous section, 
one can of course take a linear combination of each solution with the unit matrix 
to obtain another solution. 



case II: b = u®v. Here 47 takes the form 



u (g> v[v, u] 8 v = u <g> [i;,u]'it <g> v . (57) 
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A simple way to satisfy this is to set 



0. 



(58) 



So, using any two commuting matrices u and v one can construct a solution 
to (47). If the elements of one of these matrices are nonnegative, and the 



nondiagonal elements of the other are also nonnegative, then the nondiagonal 
elements of b are nonnegative. If 



(s\u — (s\v 



(59) 



then b satisfies (J9j) as well. Of course, having found a solution of this type one 
can use a linear combination of it with the unit matrix as another solution. 
Solutions 1, 4, 7, 14, 17, 20, 21, 22, 25, 26, and 28 of |§] are of this kind. 

It is possible to have other solutions to j57|). In this case, let's also use (||). 
This shows that one may rescale u and v so that ( j59"l ) is satisfied. One then 
arrives at 



u = u 



v 2 = v 



vuv — uvu, 



(60) 



if [u, v] 7^ 0. From these, it is seen that u, v, 1 — u, and uvu are projections. 
Moreover, 



(1 — u)uvu = uvu{\ — u) = 0. 



(61) 



This shows that 1— u and uvu can be simultaneously diagonalized. The diagonal 
form of them will be 

/0 N 
l-it= 
\0 1. 



'I N 


,0 0; 



(62) 



Here the elements of the above matrices are matrices themselves, and 1 is the 
unit matrix of the appropriate dimension. Writing an ansatz for v: 



"11 


"12 


"13 


"21 


"22 


"23 


"31 


"32 


"33 



(63) 
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and putting it in (pOj), one finally arrives at the following forms for u and v. 


















1 
























v° 








( V 




f 1 
















!(> 








w' 


1 





v° 








0/ 



(64) 



where all of the entries in the above matrices are matrices, and w and w' should 
satisfy 

to' = w'w = 0. (65) 

Each of the diagonal blocks of these matrices may be zero dimensional, except 
the first. It should, at least, be one dimensional. The reason is that u and v have 
at least one common left eigenvector, (s\, corresponding to the unit eigenvalue. 
Also the dimension of each block of u is equal to that of the corresponding 
block in v. Also note that if the dimension of u and v is 2 (there are two kinds 
of particles) then ther will no space left for w and w' , and u and v must be 
commuting. 

The final result is that in two dimension no new solution exists (u and v 
must be commuting), and in more than two dimensions, u and v must be of 
the form ( p4|). Of course any similarity transformation on ( |64] ) gives another 
solution to (|47|). In fact, one has to use a similarity transformation to make (s| 
a left eigenvector of u and v with unit eigenvalue. 

Two very simple subcases are b = 1 ® v and b = u <£> 1. These describe 
reactions 

A a A fi -> A a A s , with rate v 5 p , (66) 

and 

A a A^ — > A' 1 A 13 , with rate uj, (67) 

respectively. That is, in each case only one of the particles change, and the rate 
of change is independent of the type of the other particle. 
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